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APPROXIMATE  EQUATIONS  FOR  TRANSPORT  COEFFICIENTS  OF 
MULTICOMPONENT  MIXTURES  OF  NEUTRAL  GASES 
I .  Introduction 

A  description  of  multicomponent  gas  mixtures  which  contains  the 
effects  of  fluid  dynamics  and  chemical  reactions  is  necessary  for  model¬ 
ling  many  reactive  flow  systems.  A  Quantitative  description  of  such 
systems  often  must  Include  an  accurate  treatment  of  transport  phenomena, 
in  which  heat  and  particles  diffuse  through  the  mixture.  Our  particular 
Interest  in  accurately  modelling  combustion  processes  has  thus  led  us  to 
survey  and  compare  existing  theories  of  transport  phenomena  in  mixtures 
of  neutral  gases  in  order  to  Identify  representations  which  are  suitable 
for  numerical  calculations.  In  this  paper,  we  present  the  results  of  our 
work  on  (ordinary)  diffusion,  viscosity,  thermal  conductivity,  and  thermal 
diffusion. 

To  show  how  these  transport  processes  enter  our  calculations,  we 
consider  the  equations  for  conservation  of  mass,  particle  number,  momentum, 
and  energy,  which  may  be  written  as  follows: 

2_e  +  i  •  ( pv.)  -  o  (i.i) 

3 1 

-tr’-z-  - 1  •  +  \  -  Vi 
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ft  (p^0) 


-  7  •  (  pv0v0  )  -  7  •  P 


ff  -  -  l  •  (Evq)  -  7  •  (vQ  •  ?)  -  7  •  q 


(1.3) 


(1.4) 


Here  v^  is  Che  fluid  velocity  and  p, p  ,  and  E  are  the  total  mass, 
momentum, and  energy  densities  of  the  gas  mixture,  respectively.  The 
sets  of  variables  {n^},  {Q^},  and  {L^}  refer  respectively  to  the  number 
densities  of  the  chemical  species  (labelled  by  i)  and  the  rates  of 
production  and  loss  of  the  species  as  a  result  of  chemical  reactions. 

As  we  show  in  Appendix  A,  the  transport  coefficients  enter  the 
conservation  equations  through  the  pressure  tensor  P,  the  heat  flux 
3,  and  the  diffusion  velocities  (V^},  which  represent  the  average 
velocities  of  the  species  relative  to  the  fluid  velocity,  v^.  The 
pressure  tensor  P  depends  on  the  scalar  pressure  P  and  the  viscosity 
of  the  mixture.  The  variables  {V^}  and  3  are  functions  of  (ordinary) 
diffusion,  thermal  conductivity,  and  thermal  diffusion.  In  Appendix  A, 
we  consider  these  relationships  in  detail. 

The  remaining  sections  of  this  paper  present  equations  for  the 

transport  coefficients  which  represent  the  above  processes  in  the 

conservation  equations.  In  developing  and  identifying  the  formulas 

we  have  attempted  to  strike  a  balance  between  accuracy  and  computational 

simplicity.  In  the  accompanying  discussions,  we  have  supplemented  the 

1  2 

major  works  on  the  theory  of  transport  phenomena  in  gas  mixtures  ’  by 
citing  publications  which  have  extended  the  theory  and  by  commenting  on 
the  consistency  among  and  limitations  of  the  various  formulations.  We  have 
also  sought  to  clarify  points  of  confusion  arising  from  differences  in 


methods  of  approximation  and  in  notation  used  by  various  authors. 

Because  the  text  will  often  mention  references  1  and  2,  we  will 
use  the  respective  designations  "HCB"  and  "CC”  to  identify  them  for  the 
reader.  We  will  employ  the  notation  of  HCB  except  when  we  are  dealing 
with  notational  differences  among  authors  and  when  we  discuss  vectors 
(e.g.,  V),  tensors  (e.g.,  B) ,  and  the  thermal  diffusion  coefficients 
{D^} .  In  references  to  HCB  and  CC  we  cite  the  equation  and  page 
numbers,  respectively,  abbreviating  the  word  "equation(s)"  by  "eq." 

(e.g.,  HCB,eq.  (8.1-1),  516).  We  will  express  physical  quantities  in 

CGS  units  and  will  write  units  in  parentheses  beside  the  related  definitions 

and  equations. 
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II.  Basic  Formulations  of  the  Theory  of  Transport  Phenomena 


Before  stating  the  equations  for  transport  coefficients  in  multi- 
component  gas  mixtures,  we  will  outline  the  approaches  of  HCB  and  CC  in 
deriving  the  fundamental  relations  from  which  our  equations  originate. 
HCB  and  CC  begin  with  the  Boltzmann  distribution  functions  {f^}  which 
satisfy  the  Boltzmann  equations  for  the  species  in  the  gas  mixture: 


(?i 


(II. 1) 


(i  -  1,  2 . ,v). 


where  f^  is  a  function  of  the  position  vector  r  (cm) ,  the  velocity 
vector  of  the  itl1  species  v^  (cm  sec  ^) ,  and  the  time  t(sec).  The 
external  force  on  each  molecule  of  species  i  is  X^  (dynes) ;  the  mass  of 
a  molecule  of  type  i  is  m^g) ;  the  t J (f ^ >  f^)}  are  collision  integrals; 
and  v  is  the  number  of  species  comprising  the  gas  mixture.  The  authors 
assume  that  the  {f^}  may  be  expressed  in  terms  of  an  infinite  perturba¬ 
tion  series 


[0] 


+  f. 


[1] 


+  f. 


[2] 


(i  -  1,  2, . v) 

where  the  lowest  order  term  is  Maxwellian,  i.e., 


(II. 2) 


/^V/2  «p 

\  2irkT  ) 

(  2kT  J 

(II. 3) 
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In  eq.  (II. 3),  n^  is  the  number  density  of  species  i(cm  ),k  is 


Boltzmann's  constant  (■  1.380662  x  10  ^  erg  °K  ^);  T  is  the  temperature 
(°K);  and  Vg(r,t)  is  the  mass  average  (or  fluid)  velocity  of  the  gas 
mixture  given  by 


I0  (J-C>  -  p 


V 

Z 

j-1 


P3  h 


(cm  sec  1) 


(II. 4) 


In  eq.  (II. 4),  p(r,t)  is  the  total  mass  density  of  the  mixture  and  is  equal 
to  the  sum  of  the  species  mass  densities  {p^}: 


v  -3 

P  ■  Z  p  (gm  cm  )  , 

i-1 

where  the  individual  mass  densities  are  given  by 


(II. 5) 


°i  *  ni“i  (i  "  x’2 . v>- 


(II. 6) 


Also  in  eq.  (II. 4)  we  have  v^,  the  average  velocity  of  species  i,  defined 
by 


^(r.t) 


,c) 


v  dv 
-i  -i 


(cm  sec 


(II. 7) 


Given  sufficiently  rapid  convergence  of  the  series  in  eq.(II.2),  we 
may  use  the  first  two  terms  to  approximate  f^;  HCB  and  CC  use  the  following 
form: 

ff(r,Vf  ft)  *  (E.^.t)  {1-HJ^  (r,v^,t)}, 

(i  -  1,2, . v),  (II. 8) 
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where  the  function  0^  is  found  to  be  linear  in  7 taX,  the  tensor 
and  a  set  of  quantities  {d^},  so  that 

3  \> 

0  -  -  (A.  •  VfiriT)  -  l  (B, )  .  (VO-  +  n  2  C.(j)  •  d. 

i  -i  -  ot,S“l  2to  — 0  3a  -i  -j 

(i  -  1,2 . v) . 

(3) 

In  eq.  (II. 9),  and  the  set  {C^  }  are  vectors,  is  a  second-rank 


(II. 9) 


tensor,  and 


h  ■  1\ 


*  (”i  -  T>  v  -  h  -  ,1  “j  Sjl  <-'1' 


(II. 10) 


(i  -  1,2, 


where  P  is  the  total  scalar  pressure,  xi  is  the  mole  fraction  of  species  i  given 


i  n  ’ 

and  n  is  the  total  number  density  equal  to 

v 

n  «  2  n,. 

3-1  J 


(II. 11) 


(11.12) 


The  quantities  A^,  J  ,  and  B^  are  all  functions  of  the  "peculiar  velocity" 
V^,  defined  by 


-i  ^*2*^  =  -i  ~  Yq  (cm  sec  _i) 


(11.13) 
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or,  alternatively,  the  reduced  velocity. 


We  then  have 


/  Y*  -1 

u  «  ( — i.  1  V,  (cm  sec  ). 

-i  \2kT  /  -i 

(11.14) 

-i  “  W-i’ 

(11.15) 

Ct(j)  -  Ci(j)(Wi)  Wif  and 

(11.16) 

h’  »,«,)  -  i  0)  , 

(11.17) 

where  the  unit  tensor  is 


U  » 


3 

Z 

i-1 


-i-i 


10  0 
0  10 
0  0  1 


and 


(11.18) 


(11.19) 


We  must  distinguish  between  the  peculiar  velocity  ^  in  eq.  (11.13)  and 
the  "diffusion  velocity"  which  is  the  average  peculiar  velocity  of 
species  i  (i.e,  the  average  rate  of  flow  of  species  i  relative  to  the  mass 
average  velocity  of  the  mixture) .  Thus  we  have 


and,  by  eqs.(II.4)  and  (II. 6),  the  set  {V^}  satisfies  the  relation 

v  V 

E  V  -  £  p  V  *  o  (11.21) 

i-1  11-1  i=l  1  _1 

From  the  above  formalism,  HCB  and  CC  derive  the  transport 
coefficients  in  terms  of  integral  equations  containing  f^^  and 
along  with  \  or  depending  on  the  particular  transport  coefficient 

being  calculated.  The  approximate  solution  of  these  complicated  equations 
1  requires  two  further  steps: 

(1)  Expand  the  functions  A^,  and  B^  in  terms  of  the  Sonine 

polynomials  {S  (m^(x) } ,  which  satisfy  a  particularly  convenient  orthogonality 

3 

condition(see  Appendix  B) .  This  gives  us  the  form: 

t.(j)  (lO  -2  t  (J}  S™  (W  2)  (11.22) 

i  i  im  n  i 

IQ 

where  t^^  is  equal  to  A^,  B^,  or  ^  and  the  number  n  is,  respectively, 
3/2,  5/2,  or  3/2. 

(2)  Derive  (integral)  equations^’^  satisfied  by  the  sets  {A^}, 

{B^} ,  and  {C^^}  and  use  them  in  conjunction  with  (11.22)  to  solve  for 

expansion  coefficients  t . ^ .  These  steps  result  in  approximate  equations 

im 

for  the  transport  coefficients  in  terms  of  a  finite  set  of  the 
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wmmmm 


Unfortunately  for  the  reader,  HCB  and  CC  calculate  the  expansion 
coefficients  t^^  in  different  ways  .  To  develop  an  "wt-  order"  approxi¬ 
mation,  HCB  initially  use  expansions  consisting  of  only  the  m  lowest  order 
Sonine  polynomials,  solving  variationally  for  the  coefficients.  CC  use 
expansions  containing  an  infinite  set  of  Sonine  polynomials  and  develop 
expressions  which  may  in  theory  be  solved  for  the  infinite  set  of  coef¬ 
ficients.  They  then  cut  off  the  expansions  at  order  m  to  obtain  their 
"m^  -  order"  approximations  for  the  transport  coefficients.  These 
techniques  result  in  somewhat  different  equations  at  each  order  of 
approximation;  however,  we  demonstrate  later  that,  for  most  purposes,  the 
differences  are  not  significant. 

The  approximate  equations  for  the  mixture  transport  coefficients  which 
are  derived  by  HCB  and  CC  take  the  form  of  ratios  of  determinants  of  order 
£  s  v  +1  and  £-1,  respectively,  where  again  v  is  the  number  of  species  in 
the  mixture.  The  elements  of  the  determinants  are  complicated  functions 
of  collision  integrals,  the  evaluation  of  which  would  require  a  significant 
amount  of  computer  time.  Consequently  the  expressions  derived  by  HCB  and  CC 
are  impractical  for  use  in  models  of  reactive  flows.  Subsequent  work, 
however,  has  produced  accurate,  semi-empirical  equations  which  are  simple  to 
evaluate.  The  remaining  sections  of  this  paper  discuss  these  equations. 
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III.  (Ordinary)  Diffusion 


A.  Binary  Mixtures  (v=2) 

We  will  use  the  symbol  £  to  denote  the  "ordinary"  (or"concen- 
tration")  diffusion  coefficient  of  a  binary  mixture  containing 
species  i  and  j.  Marrero  and  Mason^  have  found  that,  for  many 
binary  mixtures  of  dilute  gases,  the  existing  experimental  data  are 
sufficiently  extensive,  consistent,  and  accurate  to  permit  the 
development  of  a  semi-empirical  expression  which  describes  the  variation 
of  over  a  wide  temperature  range (300°K  i  I  s  10000°K)  and  which  has 
a  standard  deviation  of  10%  or  less.  Based  on  the  qualitative 
temperature  dependence  of  the  quantity  P£^  ,  where  P  is  the  pressure, 
Marrero  and  Mason  used  the  following  form  to  fit  the  more  accurate  and 
consistent  data: 


pAemP 

ij 


AT 


(III.l) 


where  A,  s,  cp^,  S  and  S’  are  empirical  constants.  The  empirical 
diffusion  coefficient,  *-s  in  units  of  cm^  sec  \  and  P  is  in 

6  *2 

units  of  atmospheres. (One  atmosphere  equals  1.01325  x  10  dynes  cm  .) 

The  terms  containing  S  and  S'  are  the  "Sutherland-Reinganum  terms",  which 
are  used  to  account  for  the  attractive  portion  of  the  intermolecular 
potential.  Marrero  and-  Mason  fit  most  of  the  less  precise  data  to  the  form 
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(HI. 2) 


P  S 


emp 


A'T 


where  A'  and  s'  are  empirical  constants.  The  function  in  eq.(III.2) 
is  much  easier  to  evaluate  and  provides  sufficient  accuracy  for  reactive 
flow  calculations.  In  those  cases  for  which  Marrero  and  Mason  have  fit 
the  data  to  eq.(III.l),  we  suggest  using  eq.(III.l)  to  generate  values  of 
at  various  temperatures  and  then  fitting  these  to  eq.(III.2). 

For  most  intermediate  or  highly  reactive  species,  experimental  data 
do  not  exist  over  the  large  range  of  temperatures  required  for  reactive 
flow  calculations.  We  must  then  rely  on  a  theoretical  expression  for 
<0^  which  is  based  on  kinetic  theory.  HCB  and  CC  have  derived  the 
i  following  equation  for  ]^,  the  first  (or  lowest  nonzero  order) 
approximation  to  f) 


ir  _ 

yT3(M  +  M  ) 
3 

2MiMJ 


(cm2  sec-1)  (III. 3) 


where  is  the  molecular  weight  of  species  i,  P  is  the  pressure  in 
atmospheres,  an<*  often  called  the  collision  diameter,  is  a  force 

constant  in  the  potential  function 


rij 


£  f  (!it 
«  rij 


(III. 4) 


which  describes  the  interactions  between  molecules  of  species  i  and  j . 
We  may  describe  o^  as  the  value  of  the  spherical  coordinate  r^  at 
which  is  equal  to  zero.  The  other  constant.  E^y  equals 
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the  depth  of  the  potential  well(or  the  maximum  energy  of  attraction). 


Unfortunately  values  of  o^  ande„  are  not  available  in  most  cases; 
so  we  will  use  the  customary  representations. 


°1J  *2  (ol  +  V- 


(III. 5) 


-8 

measured  in  A  (1A  *  10  cm) ,  and 


(III. 6) 


measured  in  ergs.  In  eq.(III.5)  and  (III. 6),  the  symbols  a  and  e  ^ 
denote  the  force  constants  in  eq.  (III. 4)  for  interactions  between  molecules 
of  the  same  species  i.  Eq.(III.5)  is  exact  for  rigid  spherical  molecules 

g 

while  eq.(III.6)  follows  from  the  theory  of  the  "London  dispersion  forces"  . 

* 

Finally,  the  collision  integral  (T*^),  is  actuaHy  a  ratio^  ; 


0 


(1,1) 

ij 
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,(1,D 

.11 


[a 


(i,D 

ij 


«  1 

3 

Rigid  Sphere 


(III. 7) 


where  the  collision  integral  is  related  to  that  defined  in  CC  by 


Q(1,1>  -  n(1)(i). 
ij  j 


(III. 8) 


(1  1)  (1  1) 

Note  that  ’  measures  the  departure  of  *  from  its  value  in  a 

model  assuming  rigid  spherical  molecules.  The  collision  integral  is 

a  function  of  a  reduced  temperature 
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(HI. 9) 


„  *  kT 
T  *  — 


To  obtain  values  of  e  and  o„  we  suggest  using  the  tabulations  of 

and  given  by  Svehla^  and  HCB^  for  interactions  described  by  the 


Lennard-Jones  potential, 


$  (r) 
LJ 


^>12  -&] 


(111.10) 
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HCB  also  provide  a  table  of  collision  integrals  0  '  for  small 
values  of  £  and  s.  Table  IIX-1  compares  values  for  %  calculated  from 
the  semiempirical  formulas  of  Marrero  and  Mason  with  those  calculated  Irom 
eq  (III.3-1Q).  At  low  temperatures  the  agreement  is  quite  good;  however, 
the  disparity  increases  with  temperature.  As  in  the  case  of  eq.  (1II.1), 
we  may  transform  eq.(III.3)  to  the  form  of  eq.(III.2)  by  using  eq.(111.3) 
to  generate  values  of  PJD^  at  various  temperatures  and  then  fitting  these 
values  to  eq.(III.2). 


Based  on  the  assumption  of  rigid  spherical  molecules,  Cheung  et  al. 
14 

and  Mason  and  Saxena  have  derived  an  approximate  equation  for  for 
mixtures  of  nonpolar  gases: 


13 
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(III. 11) 
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In  eq.  (III. 11),  A®  (erg  ca  *  sec  1  °K  S  is  the  coefficient  of  thermal 
conductivity  calculated  for  a  gas  consisting  only  of  species  i  at  the 
same  pressure  and  temperature  as  the  mixture  being  studied  (T^  *  T^^ , 

Pi  »  Pm<^)  and  with  the  internal  degrees  of  freedom  "frozen".  The 
last  condition  means  chat  species  i  is  treated  as  a  monatomic  gas 
when  determining  A°  (Section  V).  Similarly,  (gm  cm”1  sec”1) 
is  the  coefficient  of  viscosity  of  a  gas  consisting  of  only  species 
i  at  the  same  temperature  and  pressure  as  the  mixture  being  studied 
(Section  IV).  The  quantity  is  the  coefficient  of  "self-diffusion", 

for  which  we  can  use  the  first  approximation,1^ 


LP  ] 

^iiJl 


0.0026280 

pc  2n(1,1)tT*) 
i  ii  ur 


,  2  -In 

(cm  sec  ) , 


(III. 12) 


where  all  quantities  are  defined  as  in  the  case  of  eq.  (III. 3)  and  we 
have  substituted  and  T*  for  and  T*^,  respectively.  In  eq.  (IV. 10) 
we  also  present  a  semi-empirical  expression  fox  Table  III-2 

compares  experimental  values  ofe $  with  those  calculated  from  eq. 

(III. 12).  Again  the  agreement  is  quite  good. 
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After  HCB,  Table  8.4-13,  581 


12 

Cheung  et  al,  have  shown  that  the  0^  calculated  frogt  eq. (111.11)  using 

experimentally  determined  viscosity  coefficients  agrees  quite  well  with  the 

ratio  calculated  from  eqs.  (III. 3)  and  (III. 12)  using  a  Lennard-Jones 

potential.  Discrepencies  will,  however,  arise  because  the  derivation  of  eq. 

(III. 11)  relies  on  the  assumption  of  rigid  sperical  molecules. 

B.  Multicomponent  Mixtures  ( v  >  2) 

The  calculation  of  diffusion  coefficients  (denoted  as  opposed  to 

)  lot  mixtures  of  more  than  two  species  requires  the  evaluation  of 

17  18 

determinants  of  order  v  and  the  associated  minors .  ’  For  more  than  four 

species,  the  calculation  of  the  set  {D^}  could  become  expensive;  fortunately, 
both  HCB  and  CC  have  derived  expressions  relating  ordinary  diffusion  to  other 
transport  properties  in  terms  of  the  first  approximation  to  the  binary 

diffusion  coefficient  given  in  eq.(III,3).  We  mention  that,1  while  ^  is 
the  same  in  the  formulations  of  HCB  and  CC,  the  first  approximation  to  D 


(denoted  [D^ ^)  is  not.  We  infer  this  from  the  fact  that 
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fill 


CC 
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(III. 13) 


while  eq.(III.13)  does  not  hold  in  the  case  of  HCB.^  Thus  we  have 


MCB,t  [Dij]icc 


(III. 14) 
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a  situation  which  we  also  find  for  coefficients  of  thermal  diffusion  (Section 
VI). 

From  a  practical  point  of  view,  the  above  disagreement  between  the 
theoretical  formulations  may  not  be  too  significant.  Marrero  and  Mason^ 
point  out  that  the  major  difference  between  and  lies  in  the  dependence 
of  on  the  relative  amounts  of  all  species  rather  than  just  species  i  and 
j  as  in  the  case  of  .  Because  the  experimental  uncertainties  in  measurements 
of  and  are  often  of  the  same  order  of  magnitude  as  the  variations  in 
composition,  Marrero  and  Mason  suggest  that 

emp  emp 

Dij  *  ^ij  *  (III. 15) 

where  "emp"  signifies  empirically  determined  quantities. 

C.  Summary 

From  the  above  discussions,  we  find  that  good  experimental  data  and  simple 
empirical  formulas  exist  for  binary  diffusion  coefficients  of  many  species.  For 
those  which  are  not  adequately  covered  by  measurements,  we  may  use  either  eq. 
(III.3)or  the  combination  of  eqs.  (III. 11)  and  (III. 12).  Because  HCB  and  CC  have 
constructed  the  equations  describing  diffusion  and  other  transport  phenomena 
in  terms  of  the  binary  diffusion  coefficients,  we  need  not  calculate  the 
multicomponent  coefficients  of  diffusion,  which  are  defined  differently  by 
HCB  and  CC. 
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IV.  Viscosity 
A.  Pure  Gas 


For  a  pure  gas  (i.e.,  one  which  contains  only  one  species),  we  may 

20 

calculate  the  coefficient  of  viscosity,  r\^,  by  using  the  expression 


266.93  x  10 
2  (2,2)  *.  * 
ai  nii  (V 


r7  /v 


,  -1  -1. 
(g  cm  sec  ) 


(IV.l) 


where  ^denotes  the  first  (i,e.,  lowest  nonzero  order)  approximation  to 

(2  2^  9 

n^.  In  eq.  (IV.l),  *  is  a  ratio  of  collision  integrals  defined  by 
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(2,2)*.  _^i 


(2,2) 


ii 


[aj2>2)] 


«  1 


(IV. 2) 


Rigid -Sphere 


(2  2) 

where  the  collision  Integral  0^  *  is  related  to  that  defined  in  CC  by 


A  ,  ^(2)  ^) . 

ii  'ii  K  ) 


(IV. 3) 


HCB  have  compared  experimental  viscosity  data  with  values  calculated  from 
eq. (IV.l)  and  eq.  (III. 10)  for  several  species  over  a  temperature  range  of 
80  -  1500°K  and  have  found  deviations  of  less  that  5%  in  all  cases.  For 

sufficiently  restricted  "temperature  ranges  we  might  simplify  the  calculations  by 

(2  2)  * 

substituting  a  constant  for  0^  *  ,  which  is  a  slowly  varying  function  of 

reduced  temperature  T  .  As  we  showed  in  eq.  (III.  11),  [n^]^  and  are 

closely  related;  we  will  show  this  relationship  explicitly  in  Section  V. 
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B.  Multicomponent  Mixtures  (v  *  2) 

For  mixtures  of  two  or  more  gases,  we  may  express  the  first  approximation 

to  the  coefficient  of  viscosity,  Ch^ x ,  as  a  ratio  of  determinants  of  order 

v  +  1  and  v,  respectively,  in  which  the  elements  are  complicated  functions  of 

22 

collision  integrals.  Of  more  use  in  reactive  flow  calculations  is  an 

alternate  expression  derived  by  expanding  the  determinants  and  discarding  the 

23 

off-diagonal  elements.  This  results  in  the  functional  form 
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i»l 


k*i 


(IV. 4) 
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l  +  o 


Chill  RT 
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k-1 
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*k 


^ik’l 


where  Chili  is  the  first  approximation  to  the  coefficient  of  viscosity  of 

a  gas  containing  only  species  i  at  the  temperature  and  pressure  of  the 

mixture  being  studied,  o  is  an  empirical  constant,  and  R  is  the  "molar  gas 

o 

constant"  given  by 


R  -  Nk  -  8.31441  x  10  7  erg  mol'1  °K_1 


(IV. 5) 


with  Avogadro's  number  given  by 


6.022045 


23  -1 

X  10  mol  1 


(IV.  6) 
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We  oay  also  express  eq.  (IV. 4)  in  terms  of  P^,  the  density  (g  cm  )  of  a 
gas  containing  only  species  i  at  the  temperature  and  pressure  of  the  mixture, 
by  using  the  ideal  gas  law: 


PV  -  TtRT.  (IV. 7) 

Here  is  the  number  of  moles  of  a  gas  containing  only  species  i  at  the  same 
temperature  and  pressure  as  the  mixture  (and  not  the  number  of  moles  of  species 
i  in  the  mixture).  We  then  obtain 
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i-l 


1  + 


°i  *i 
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k-i 

k*i 


*k 


^iki 


(IV. 8) 


We  point  out  again  that  p'  is  not  the  density  of  species  i  in  the  mixture; 

rather  p'  is  Cue  density  corresponding  to  T,  in  eq.  (IV. 7).  Using  a 
i  1 

24 

heuristic  treatment,  Buddenberg  and  Wilke  obtained  eq.  (IV. 8)  with 

H  .  and  n.,  substituted  for  [n  .  ],  and  [n4],  and  evaluated  a  to  be  approximately 
mix  i  mix  1  11  o 

1.385  by  comparing  (IV. 8)  with  available  viscosity  data  for  binary  mix¬ 
tures.  Table  IV-1  shows  that  we  can  calculate  values  of  n  which  are 

mix 

in  excellent  agreement  with  experimental  data  for  mixtures  of  two,  three, 

and  four  components  by  using  eq.  (IV. 8)  with  otQ  *  1.385  and  experimental 

values  of  p'^  and  n^.  The  dependence  of  eq.  (IV. 8)  on  binary  diffusion 

coefficients,  for  which  data  are  often  unavailable,  can  make  calculations 

25 

of  mixture  viscosities  more  difficult  and  less  accurate.  Wilke  ,  there¬ 
fore,  developed  the  following  equation  for  [ v ^  in  terms  of  the  quan¬ 
tities  {0^^} ,  defined  by  eq.  (III. 11): 
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Table  IV-1 


Comparison  of  Theoretical  and  Experimental  Values  of 
Viscosity  Coefficients  of  Gas  Mixtures  at  298°K  and  1  Atmosphere 


Gas  ’ 

Mixture 

Mole 

Fractions 

1.  (Experimental) 26 
mlx10"®a;  cm~l  sec"l 

W(Tirr,ti 

10  cm 

55 

NJ 

1 

o 

N> 

0.1864-0.3136 

0.7822-0.2178 

200.8 

184.3 

197.4  24, 
180.6  24’ 

h2-C° 

0.1927-0.8073 

0.6947-0.3053 

171.7 

144.9 

171.9  24, 
149.6  24, 

He-Ar 

0.3405-0.6595 

0.7565-0.2435 

227.8 

227.0 

227.0  24> 
227.8  24’ 

Ne-H2-C02 

0.333-0.333-0.333 

185.7 

195.3  24 

Ne-H2-C02 

’CC12F2 

0.25-0.25-0.25-0.25 

168.1 

162.3  24,S 

i/. 


*  After  Buddenberg  and  Wilke. 

a.  From  eq.  (IV. 8)  of  this  paper 

b.  From  eq.  (IV. 9)  of  this  paper 
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Depending  on  the  species  in  the  mixture  and  the  size  of  the  reduced  temp¬ 
erature  vs 
replacing 


&  A 

erature  range  of  interest  (T^,  Tg) ,  we  might  simplify  eq.  17.  12  by 
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'  a, 2>* 
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(V 
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n  (t  ) 

ii  1  r 


(IV. 13) 


with  its  average  over  (T^,  Tp .  For  models  involving  a  large  range  of 

reduced  temperatures,  accurate  calculations  will  require  the  use  of  a 
(2.  * 

Table  of  ’  values  vs.  T*  or  an  accurate  fit  of  a  convenient  func¬ 
tional  form  to  Che  tabulated  values. 


V.  Thermal  Conductivity 

The  expression  for  the  coefficient  of  thermal  conductivity  for  a 

mixture  of  reacting  gases,  X  is  often  quite  complex,  as  we  indicate  in 

Appendices  A,  D,  and  E.  Depending  on  the  speeds  of  the  reactions  and  on 

the  lifetimes  of  the  excited  states  of  each  molecule,  we  may  express 

X  .  as  follows: 
mix 


mix 


X°  + 

mix  mix  mix 


(V.l) 


In  eq.(V.l),  (sometimes  called  the  "translational  thermal 

conductivity")  is  a  function  of  the  set  {A°  J  i*l,2,...,v}  of  thermal 
conductivities  of  the  individual  species,  each  treated  as  monatomic. 

The  second  term.  A*”*,  accounts  for  the  transport  of  internal  energy 
by  molecules  in  different  quantum  states,  and  the  last  term, 
describes  heat  conduction  resulting  from  changes  in  the  composition  of 
the  mixture  through  chemical  reactions.  We  must  point  out,  however,  that 
the  heat  flux  resulting  from  changes  in  chemical  composition  of  the  mixture 
may  be  represented  as 

comp  „  .  xcomp  7T 
-  mix 

only  if  the  steady  state  chemical  composition  is  very  nearly  in  equilibrium 

with  the  local  temperature  T.  Because  this  depends  on  the  reaction  rates 

in  the  forward  and  reverse  directions,  the  term  \c°mp  in  eq.(V.l-2)  is 

mix 

not  always  applicable.  We  have,  therefore,  derived  the  following  general 
expression  for  a  comP  in  Appendix  D: 
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i-1  1  -1  1 

where  71  is  the  number  of  moles  in  the  gas  mixture  and  <H^>  is  the 
experimentally  determined  enthalpy  per  mole  of  gas  species  i.  We  may 
now  restate  eq.(V.l)  as 


mix 


mix 


.  int 
mix 


(V.4) 


In  the  remaining  sections,  we  will  discuss 


\  .  and 
mix 


in  detail. 

mix 


A.  Pure  Gas 

For  a  gas  consisting  of  only  one  species,  we  first  consider 

which  is  calculated  assuming  that  the  gas  is  monatomic  or  that  the 

27 

internal  degrees  of  freedom  are  frozen.  HCB  give  the  following  expression 
for  the  first  approximation,  [1°]^  : 


t^l 


8.322  X  10 
2  (2,2)*  * 
°i  nu 


,  T  v!S  /  “1  -1  0..-1, 

(  rr-  )  (erg  cm  sec  K  ) . 

Mi 
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By  eq. (IV. 1)  and  eq.(V.5),  we  find  that  [n^ ] ^  and  [a®]^  are  clsoely  related,  i.e 


-9 


(nil1  -  3.208  X  10 


(V.6) 


For  polyatomic  molecules,  of  course,  the  internal  degrees  of  freedom 

are  not  frozen,  and  the  diffusion  of  molecules  of  species  i  in  different 

quantum  states  will  transport  energy  through  the  gas.  We  may  represent  this 

27 

process  approximately  by  using  the  Eucken  factor,  E^,  with  eq.(V.5)  as 
follows  : 


26 
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Hirschef elder  has  derived  the  following  equation  for  E^  from  the 

kinetic  theory  of  gases  (Appendix  E) : 


Ei  -  0.115  +  0.354  ^ 


-  0.115  +  0.354 


(V.8) 


where  c  .  is  the  constant  pressure  molar  specific  heat  for  molecular 
species  i. 


f2i 

Cvi 


(V.9) 


and  c  is  the  constant  volume  molar  specific  heat  for  molecular  species 
i. 


Hirschf elder  assumed  the  following  in  order  to  derive  eq.(V.8): 

(1)  The  gas  exists  in  a  steady  state. 

(2)  The  composition  of  the  gas  (i.e.f  population  of  the  various  quantum 
levels)  is  in  equilibrium  with  the  local  temperature.  Only  when  this 
assumption  holds  will  the  heat  conductivity  be  independent  of  boundary 
conditions,  permitting  us  to  express  the  heat  flux  due  to  internal 
degrees  of  freedom  as 


int  . int 

3i  "  -  xi 


7T  . 


(V.10) 


If,  for  example,  an  appreciable  fraction  of  the  molecules  were  to  exist 


27 


in  metastable  states,  the  population  of  states  would  not  be  a  function 
temperature  only.  Eq.  (V.10)  would  then  be  incorrect  and  X*nt  would  not 
provide  a  meaningful  representation  of  polyatomic  effects. 

(3)  The  coefficients  of  diffusion  for  all  quantum  states  with  a  nonnegli- 

gible  population  are  approximately  the  same.  This  is  usually  the  case, 

except  when  excited  electronic  states  are  appreciably  populated.  Such  states 

have  collision  diameters  a?  which  are  on  the  order  of  3c? ,  where  a  is 

i  i 

the  label  of  the  quantum  state  and  a  »  0  signifies  the  ground  state.  The  large 

values  of  oa  for  excited  electronic  states  result  in  much  smaller  diffusion 

CL  -2 

coefficients,  which  vary  as  (c  )  ,  than  for  the  ground  electronic  state, 

and  consequently  the  Eucken  factor  would  predict  a  X  which  is  larger  than 
the  experimental  data  would  indicate. 

(4)  Nonadiabatic  collisions  do  not  appreciably  distort  the  molecular 
distribution  functions  {f®  (r,  v^,  t)  |  a  ■  0,1,2,....}  for  the  various 
quantum  species  a.  This  holds  in  the  case  of  electronic  and  vibrational 
states.  However,  at  "low"  temperatures,  where  rotational  transitions  become 
important  in  collislonal  energy  transfer,  rather  large  distortions  do  occur 
in  the  distribution  functions,  leading  to  anomalously  small  experimental 
values  of  X^.  For  nonpolar  molecules,  "low  temperature"  means  ~170°K, 
while  for  polar  molecules,  "low  temperature"  signifies  room  temperature. 
Fortunately  the  Eucken  factor  works  well  at  higher  temperatures  which 

are  of  interest  in  reactive  flow  problems,  because  the  distortions  due  to 
rotational  transitions  are  less  Important. 

Appendix  E  shows  how  these  assumptions  are  used  to  derive  E  .  Subsequent 

32-34  i 

work  has  accounted  for  the  effects  of  the  inelastic  collisions 

mentioned  above  and  has  resulted  in  an  Eucken  type  factor  which  involves 


only  measurable  quantities  plus  the  relaxation  time  associated  with  the 

35 

transfer  of  energy  due  rotational  transitions  in  collisions.  CC 

indicate  that  these  results  are  highly  model  dependent  and  are  currently 

impractical  for  application  to  reactive  flow  calculations. 

B.  Multicomponent  Mixtures  (v  £  2) 

36  o 

HCB  derive  [X  ,  ].,  the  first  approximation  to  X  .  for  monatomic 
mix  1  rr  mix 

gases,  in  terms  of  a  ratio  of  determinants  of  orders  2v  +  1  and  2v, 

respectively,  plus  a  term  involving  thermal  diffusion  and  ordinary  diffusion 

37-41 

coefficients.  Subsequent  work  has  led  to  a  simpler  equation 


involving  only  a  ratio  of  determinants  of  order  v  +  1  and  v,  respectively. 


As  we  indicated  in  Section  II,  the  "first  approximation",  [X°n|<v]1 ,  differs 

from  author  to  author  because  different  methods  of  approximation  have  been 
37  42 

used.  Mason  used  Kihara's  method,  in  which  derivatives  with  respect  to 

f  0  a\ 

temperature  of  a  particular  set  of  functions  of  collision  integrals  ’ 

39 

are  discarded.  Muckenfuss  and  Curtiss  used  the  variational  equations  of  HCB, 
41 

and  CC  used  truncation  of  an  infinite  series  expansion  to  develop  a  set 


of  v  coupled  equations  for 


where 


[x°.  L  =  z 

1  i-i 


W)]1- 


(v.il) 


We  show  in  Appendix  C  that  the  equations  of  CC  lead  to  a  ratio  of 

determinants  of  orders  v  +  1  and  v,  respectively,  with  elements  that  differ 

37  39 

only  slightly  from  those  of  Mason  and  Muckenfuss  and  Curtiss. 


Unfortunately  the  above  formulas  depend  on  functions  of  collision  Integrals 
which  are  nonlinear  functions  of  temperature. 


14  43 

Subsequently  Mason  et  al.  ’  have  found  that  Che  off-diagonal 

elements  of  the  determinants  could  be  neglected  and  that  the  portion  of  the 

diagonal  elements  which  contains  the  collision  integrals  could  be  replaced 

25 

by  an  empirical  constant.  The  derivation  then  parallels  Wilke's  treatment 
of  the  mixture  coefficient  of  viscosity  and  results  in  the  following 
semi-empirical  formula: 
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where 


1.065  0.k 


(V.13) 


The  factor  1.065  is  the  value  of  the  empirical  constant  mentioned  above 

o 

and  is  the  only  difference  in  functional  form  between  [l^^j  and 
[ nmiac ^  1  in  e<**  Ta^e  V-l  shows  that,  for  mixtures  of  rare  gases 

(Ne-Ar-Kr  and  He-Ar-Xe) ,  we  can  obtain  excellent  agreement  with 

experimental  data  by  using  eq. (V.12 ,13)  in  conjunction  with  empirical  values 

o  14 

of  X^.  Similar  agreement  exists  for  binary  mixtures. 

We  must  now  consider  the  correction  X*-*^  due  to  the  internal  degrees  of 

mix 

30 

freedom  of  polyatomic  molecules.  Hirschfelder  has  derived  the  following 
equation: 


XiQt 

mix 
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Comparison  of  Theoretical  and  Experimental  Values  of  Thermal  Conductivity 
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After  Mason  and  Saxena,  Reference  14. 

References  for  experimental  data  are  formed  in  parentheses 

From  eq.(V-13,16)  of  this  paper. 

From  eq.(V-12,13)  of  this  paper. 


VI.  Thermal  Diffusion 

The  equations  for  the  transport  coefficients  related  to  thermal 
diffusion  are  more  complex  and  confusing  than  those  discussed  in  previous 
sections  for  the  following  reasons: 

(1)  Thermal  diffusion  is  a  "second-order"  effect,  and  to  evaluate  the 
coefficient  of  thermal  diffusion  for  species  1,  denoted  D,^,  accurately, 
we  require  a  Sonine  polynomial  expansion  with  more  terms  than  for  other 
transport  coefficients. 

(2)  The  coefficients  of  thermal  diffusion  are  more  sensitive  to  the 
composition  of  a  gas  mixture  than  are  the  coefficients  representing 
other  transport  phenomena. 

(3)  Because  of  (1)  and  (2)  above,  experiments  and  theoretical  studies 
have  often  emphasized  binary  mixtures,  thus  impairing  our  ability  to 
verify  the  accuracy  of  formulas  for  multicomponent  mixtures  in  general. 

(4)  Three  quantities  often  appear  in  theoretical  work — D^,  kTi 

(thermal  diffusion  ratio),  and  (thermal  diffusion  factor).  Authors 

differ  on  the  definition  of  D Ti  and  sometimes  (e.g.,  HCB)  fail  to  show 

explicitly  how  the  three  are  related  for  multicomponent  mixtures  with 

v  >  2.  CC  are  also  inconsistent  in  their  definitions  of  k„.  for  v  *  2 

Ti 

and  v  >  2. 

(5)  Some  of  the  most  authoritative  papers  on  the  subject  also  involve 
the  most  complex  mathematics  and  notation,  requiring  a  significant 
effort  for  most  readers  to  gain  a  working  knowledge  of  the  equations. 

In  this  section  we  deal  with  the  above  problems  in  addition  to 
providing  a  set  of  useful  equations  for  thermal  diffusion  coefficients 
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of  multicomponent  mixtures.  While  we  use  and  compare  information 
derived  from  several  sources,  we  prefer  CC  because  the  relationships 
among  {D^,  ^ }  are  stated  most  clearly  and  completely.  For  this 

reason,  we  use  the  notation  of  CC  in  this  section,  rather  than  the 
notation  of  HCB. 

A.  Thermal  Diffusion  Coefficients  {D,^} 

The  product  of  Vfa  T  and  D_. ,  the  thermal  diffusion  coefficient 
for  species  i  in  a  multicomponent  mixture,  appears  in  eq.  (A. 3)  of 
Appendix  A  for  the  diffusion  velocity  V^  Upon  reviewing  the  litera¬ 
ture,55-58  we  find  that  CC,  HCB,  and  Waldman  (W)  all  define  DTi 
differently.  Because  we  will  be  more  interested  in  the  thermal 

diffusion  ratio,  ,  which  is  identical  in  all  three  formulations, 

HCB  W 

this  subsection  will  only  explore  the  relationships  among  DTi  ,  D^.^, 

CC 

and  Dt±  rather  than  giving  expressions  for  evaluating  them  in  terms 
of  properties  of  a  given  gas  mixture.  In  the  next  subsection,  we  will 
show  how  they  are  related  to  the  set 

The  equations  defining  the  set  tD^}  in  the  three  formulations 
are  as  follows: 


1 


_  CC 

CC.  4  -  -  ^  iyij  -  Bjl 


(VI. la) 


HCB:  V, 


4  I  «  D^Bd. 

nip  j-1  j  ^ 


“i  11 


(VI. lb) 


W:  Vj 


(VI. c) 
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Eq.  (Vl.la-c)  show  explicitly  Che  relationships  of  che  respective 


to  the  diffusion  velocity  and  multicomponent  diffusion  coefficients. 

cc 

We  readily  see  that  the  associated  units  are  cm2  sec"l(D^),  g  cm-1 
sec” 1  (D^®) ,  and  cm2  sec'^D^.).  The  first  term  of  eq.  (VI. lb)  differs 


Ti 


Ti 


in  sign  from  chose  of  eq.  (Vl.la.c)  because  the  related  term  in 

0 .(£,  t)  as  defined  by  HCB  (see  eq.  (II. 9))  differs  in  sign  from 

the  corresponding  term  in  the  other  formulations.  The  quantity  is 

CC 

a  "generalized  diffusion  coefficient"  and  is  not  equal  to  ,  the 
multicomponent  diffusion  coefficient  in  the  formulation  of  CC,as  we 
show  in  Appendix  F.  Waldmann  defines  his  diffusion  velocity  V£,  as 
the  velocity  of  species  i  relative  to  the  average  particle  velocity  of 
the  mixture. 


v' 
— o 


=  v 


XirL’ 


(VI. 2) 


rather  than  the  mass  average  velocity  v  ;  so  we  have 


(VI. 3) 


This  is  the  primary  distinction  between  the  formulations  of  CC  and 
Waldmann . 5  9 

For  consistency  with  eq.  (11.21),  which  constrains  the  diffusion 
velocities  in  the  formulations  of  CC  and  HCB  ,  and  with  the  equation 

v  _ 

Z  x.  v;  -  o  (VI. 4) 

i-i  i-i 
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in  the  treatment  of  Waldmann,  the  respective  authors  have  defined 
constraint  equations  for  the  sets  {D^}: 


cc:  2  p,  Drr  -  0 

i*l  11 


HCB:  Z  =  0 


W:  Z  x.  d!,,  ■  0 
-t-i  i  Ti 


(VI. 5a) 


(VI. 5b) 


(VI. 5c) 


From  eq.  (VI. 1)  and  (VI. 5),  we  infer  that  the  values  of  given  by 
the  three  formulations  are  related  by 


„  nCC  .HCB  -W 

pi  °Ti  *  °Ti  *  xi  °Ti 


(VI. 6) 


This  will  be  of  use  to  us  in  Section  VI. B  and  Appendix  F.  For  complete¬ 
ness,  we  also  list  the  constraints  on  the  multicomponent  diffusion 
coefficients: 


CC  i  Pj  A 


i  “ij 


PiAij 


(VI. 7a) 


-HCB  - 
HCB:  Dii  -  0 


(VI. 7b) 


W:  Z  x  DW 

J  i] 


m  Z  x.  D  »  0. 

-i.i  1 


(VI. 7c) 


?rom  eq.  (VI.7a-c),  we  see  that  the  ordinary  diffusion  coefficients 
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for  multicomponent  mixtures  are  quite  different.  Fortunately,  as  we 
have  pointed  out  in  Section  III,  HCB  and  CC  have  derived  expressions 
for  the  transport  coefficients  and  diffusion  velocities  in  terms  of 
the  binary  diffusion  coefficients  which  are  equal  in  the  two 

treatments . 

A  general  equation  for  [D^l^.the  first  approximation  to  the 
thermal  diffusion  coefficient  for  species  i,  appears  in  HCB.60  In  the 
next  section,  we  will  show  how  D,^  may  be  obtained  from  the  thermal 
diffusion  ratios 


B.  Thermal  Diffusion  Ratios  {k^}  and  Thermal  Diffusion  Factors  {o  } 
Most  data  from  experimental  studies  of  thermal  diffusion  are 


presented  in  terms  of  the  thermal  diffusion  ratios  {k^},  which  are 
unitless.  We  may  understand  this  emphasis  on  kTi  by  considering  eq. 
(A. 4)  of  Appendix  A, 
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kTi  -*71  T  —  “ 
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's' 

t—i 

3-1 


x.x. 


(V, 


V- 


(A.  4) 


Usually  measurements  take  place  after  the  gas  mixture  has  reached 
equilibrium  and  under  conditions  in  which  pressure  gradients  and 
external  forces  are  negligible.  Eq.  (A. 4)  then  becomes 


-xl*  kTl-*  I  *  °> 


(V 1.8) 


and  k^  therefore  indicates  the  changes  in  equilibrium  concentration 
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of  species  i  caused  by  a  temperature  gradient  in  the  sample.  Because 
is  so  closely  related  to  experimental  measurements,  we  expect  that, 
unlike  the  situation  for  the  thermal  diffusion  coefficients  {D^},  the 
treatments  of  CC,  HCB,  and  Waldmann  should  provide  values  of  kp^ 
which  are  very  close  or  Identical.  As  an  example,  we  show  in 
Appendix  F  that  kp^  is  identical  in  the  treatments  of  binary  mixtures 
given  by  CC  and  HCB. 

To  define  the  thermal  diffusion  ratios  formally,  both  CC61  and 
Waldmann57*59  use  the  same  type  of  equation,  i.e., 


and 


Aij  kTJ 


Dij  kTJ  * 


(VI. 9) 


(VI. 10) 


respectively.  Notice  that  by  using  eq.  (VI. 6)  with  either  eq.  (VI. 9) 
or  eq.  (VI. 10),  we  may  calculate  the  thermal  diffusion  coefficients 
in  any  of  the  formulations  given  values  for  (kp^}  and  a  set  of 
diffusion  coefficients.  Because  the  thermal  diffusion  coefficients 
are  related  by  eq.  (VI. 5),  we  require  the  auxilliary  condition,92 


S  L'O  ,  (VI.  11) 
i-1 

to  define  the  set  {kp^}  unambiguously.  The  thermal  diffusion  factors 
{a^j}  ere  defined  in  terms  of  the  thermal  diffusion  ratios  by  the  set 
of  equations57,53,63 
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(i-l.  2,  •••,  v) 


(VI. 12) 


“u  ■  R  Vj“i 


with  the  auxilliary  conditions 


-aJi 


(VI. 13) 


£q41,64  have  developed  the  following  useful  equation  for  the 
first  approximation  to  the  thermal  diffusion  ratio  of  species  i, 
[l^Ti  1 1 »  which  is  based  on  the  assumption  that  the  gases  in  the  mix¬ 
ture  are  monatomic  (or  that  the  Internal  degrees  of  freedom  are 
frozen) : 


V  lyy,.,  -  xjVl) 


(VI. 14) 


The  quantity  is  the  ratio  of  collision  integrals  appearing  in 
eq.  (C.10)  of  Appendix  C  and  a^  is  the  contribution  of  the  ith 
species  to  the  mixture  thermal  conductivity  In  Appendix  C, 

we  show  that 


^mixU)1!" - 


i+i  s*. 

i  k-l 


(C.l-2) 


CC  argue  that  the  internal  degrees  of  freedom  of  the  gas  molecules 
have  a  smaller  effect  on  k^  than  on  .  Monchick,  et  al.f5  add 
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the  condition  that  the  mass  and  size  differences  among  the  species 
should  not  both  be  small.  Otherwise,  inelastic  effects  due  to  non- 
spherical  interactions  (e.g.,  rotational  relaxation)  must  be  taken 
into  account.  An  example  of  a  mixture  for  which  internal  degrees  of 
freedom  may  not  be  ignored  is  the  D^-HT  system,  an  "isotopic"  mixture, 
in  which  both  molecular  species  are  approximately  the  same  in  mass 
and  size.  For  such  mixtures.  Monchick,  et  al.,  have  derived  the  correc¬ 
tions  to  due  to  the  internal  degrees  of  freedom  of  the  molecules. 

We  may  convert  eq.  (VI. 14)  to  a  more  useful  form  by  using 


1  +  C 


ij 


6. 

5 


(C.10) 


the  ideal  gas  law. 


P  -  nkl. 


(VI. 15) 


and  the  relationship  between  {n^}  and  {M^}  , 


M. 


m^  +  mj  +  Mj 


(i,  j  -  1,  2,  v), 


(C.9) 


We  then  obtain 


v  ,  xJx.MJM 

IkTi]l  *  ^  nkUflr^  (5  Cij  "  30 

v  u,.(6C*  -  5) 


1*1  W  aj  _  ai  \ 
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J-1VJ  Skn^l. 


(VI. 16) 


with 


VI 


ij 


M.M. 
1  J 


Mi  + 


M, 


(VI. 17) 


This  is  the  same  form  of  equation  derived  by  Monchick,  et  al. ,66 

for  however,  their  equation  for  includes  effects  related  to 

internal  degrees  of  freedom  and  is  more  complicated  from  a  computational 

standpoint.  Monchick,  et^  al . ,  point  out  that  their  result  reduces  to 

that  of  CC  for  mixtures  of  monatomic  gases.  The  presence  of  ,  which 

is  of  order  unity,  increases  the  complexity  of  eq.  (VI. 16).  If  the 

*  * 

reduced  temperature  range  of  interest  (T  ,  T  )  is  small  enough,  we 
might  simplify  our  calculations  by  substituting  the  average  value  of 
over  (T*,  T*)  for  (T*)  in  eq.(VI.16)  : 

* 

T 

*  1  ,  B  *  *  * 

Ctj  -  t4  /*  VT  )  dT  •  ^VI-18> 

tb  “  ta  ta 

For  large  temperature  variations,  accurate  calculations  will  require 

*  * 

the  use  of  a  table  of  values  vs.  T  or  an  accurate  fit  of  a 
convenient  functional  form  to  the  tabulated  values. 

While  eq.  (VI. 15-16)  should  give  an  adequate  representation  of 
thermal  diffusion  for  reactive  flow  calculations,  we  should  point  out*7 
that  past  calculations  of  for  binary  mixtures  have  not  agreed  as 
closely  with  experimental  values  as  have  similar  calculations  of  the 
other  transport  coefficients.  Thus  the  reader  may  require  higher-order 
approximations  for  other  applications. 
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C.  Additional  Remarks 


Our  study  has  shown  that  current  literature  on  the  theory  of  thermal 
diffusion  is  not  uniform  in  the  definition  of  thermal  diffusion  coeffi¬ 
cients  and  that  treatments  of  multicomponent  mixtures  are  quite  complex. 
To  produce  the  above  discussion  of  the  topic,  we  have,  therefore, 
often  been  forced  to  put  together  fragmentary  information  coming  from 
several  sources,  and  our  arguments  have  often  been  heuristic  rather  than 
rigorous.  Where  possible,  we  have  tested  the  consistency  of  the  various 
formulations,  and  as  Indicated  in  Sections  A  and  B  and  Appendix  F,  we 
have  found  them  to  be  in  agreement. 

Two  theoretical  treatments  which  we  have  not  mentioned  are  mean- 
free-path  theories68-72  and  the  phenomenological  approach  of  van  de  Ree, 
et  al. ,  based  on  the  thermodynamics  of  Irreversible  processes.73-4 
The  mean-free-path  theories  are  only  qualitatively  useful  because  of 
the  difficulty  in  computing  accurate  mean-free-paths  related  to 
"number  density  transfer"  and  "mean  thermal  speed  transfer.”  The 
formulas  derived  by  van  de  Ree,  et^  al . ,  are  equivalent  to  the  treatment 
of  Monchlck,  et  al. ,6 5,66 in  the  case  of  binary  mixtures;  these  equations 
might  thus  be  of  use  in  detailed  studies  of  thermal  diffusion. 
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VII.  Summary  of  Equations 


For  ease  of  reference  we  summarize  the  theoretical  equations  of 
major  interest  In  evaluating  transport  properties  of  multicomponent 
mixtures  of  neutral  gases.  The  definitions  of  the  symbols  appear  in 
the  body  of  the  report  on  the  pages  identified  in  Section  VIII  "List 
of'Symbols". 


A.  (Ordinary)  Diffusion 
1.  Self -Diffusion 


r  n  ,  ,  0.0026280  / T3 

1  *  Mi 


t  2  “In 
(cm  sec  )  , 


(III. 12) 


2.  Binary  Diffusion 


*uJi 


0.0026280 

Pctij  « 


/t3(M,  +  M  ) 

80 _  /  1  3 


,  2  -1, 
(cm  sec  ) 


(III. 3) 


B.  Viscosity 

1.  Pure  Gas 


'"iii 


266.93 

o2  Q  ^2 

i  ii 


x  10~7  ^/m.T 

’2)*(T?)  > 


(g  csT^sec-1) 


(IV. 1) 
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VIII.  List  of  Symbols 


The  following  is  an  alphabetic  list  of  symbols  with  the  number  of 
the  page  on  which  each  symbol  is  defined: 

(a)  A^,  6  A^,  7  A,  10  A,  11  20  aij*  ^ 

a  ,  39  a,  28  c^,  58  A  ,  63  A*^  65 


(b)  B  i,  6 


V  7 


V  63 


BiV  65 


(c)  C1(j),  6 


7  CC,  3  cpi,  27  cvl,  27 


Cij’  39  ciV40 


(d)  V,  1  d^,  6  »  3-0  ^ij^  1’  11  0^ii,  13 

^iiV 15  V 17  [Dij]i’ 17  -« 


D.r'p,  is 


dt±,  34  D  ,  36 


Ti 


75 


6fi’  81 


V  34 


(e)  E,  2  ^  ^i’  22  Oj*  ^  i’  3? 


W  23  V  27 


e,  58 


(f)  Et.  4 


,[°1 


f'11.  4 


(g)  Tt.  27  Slk,  30  0“,  59  It®.  67  ,  70  Jt>  71 


«1.  78 


20 
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wmmmum 


(h)  HCB,  3  h£,  58  <H^  ,59 

(j)  J,  4 

00  k>  *  ^ri»  37  tkT1]1,  39  i^mix *  57 

(l)  h,  1  x|,  13  Xnll,  25  1°^,  25  x‘£,  25 

CT-  25  tXjlj,  26  Xt,  27  xfe,  27  [X^lj.  29 
£Cx(i)]l’29  Lij»  63  lL'1 , 63  Ltj,  64  I L |  ,  64 

(m)  4  M^,  11  y^,  41 

(n)  ni,  1  v,  4  n,  6  8,  20  71 ,  26  n“,68 

(o)  n^’15  ,  u  i2  a^u),  12  n£2,2)  .  19 

19  ^2)(2),  19 

(p)  &,  2  P,  6  01,6  *Q.  10  ♦y.U  13 

«Lk*  13  P»ix,  15 


(q)  Qj  »  1 


a*  2 


2COmp,  25 
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Appendix  A.  Transport  Properties  and  the  Conservation  Equations 


Below  we  summarize  and  discuss  the  equations  for  conservation  of  mass, 
particle  number,  momentum,  and  energy  of  multicomponent  gas  mixtures  under 
the  assumption  that  external  forces  and  radiation  processes  are  neglibible. 
These  equations  form  the  basis  for  models  of  reactive  flows,  and  we 
shall  see  that  the  transport  phenomena  enter  the  calculations  through  the 
conservation  equations  for  particle  number,  momentum  and  energy.  For 
detailed  derivations,  we  refer  the  reader  to  the  treatments  of  Williams^ 
and  Landau  and  Lifshitz.^ 

The  "continuity  equation"  which  represents  conservation  of  total  mass, 

|6-  +  7  •  (pv  )  *  0  (A.l) 

The  transport  coefficients  do  not  appear  in  Eq.  (A.l)  because  the  related 
transport  phenomena  do  not  directly  affect  the  conservation  of  total  mass. 
However,  transport  phenomena  will  affect  the  conservation  equations  for 
mass  and  particle  number  of  each  species  i  treated  separately ,  since  the 
velocity  ^  may  differ  from  the  fluid  velocity  v q.  We  may  express  the 
conservation  of  particle  number  for  species  i  as 

1  •  -  Vl 

(A.  2) 

V  •  (n^)  -  7  •  (n^)  +  Q±  -  L^, 


is 
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where  and  refer  to  chemical  production  and  loss  processes  for  species  i. 

In  the  formulation  of  HCB,^  the  diffusion  velocities  {V.}  are  directly 

— i 

HCB 

related  to  the  coefficients  of  ordinary  diffusion,  {d“  },  and  thermal 

diffusion,  {D^B},by  the  equation 


-  n2  -HCB  .  1 _  -HCB  „  „  _ 

V.  *  -  23  m.  D.  .  d. - - —  D_.  V  2m  T. 

— i  n,p  .  ,  i  ii  —i  n.m.  Ti  — 

i  j-1  J  J  J  l  i 


(A.3) 


We  have  used  the  superscript  "HCB"  because  other  authors  define  these 

coefficients  differently;  Section  VI  shows  the  equations  of  CC  and  Waldmann 

which  correspond  to  (A.3).  We  may  also  calculate  {V^}  from  the  set  of  v-1 

78 

independent  equations 


<L  +  k  V  in  T  -  -L  (£,  "  I  ) 

j-1  d±y  ^  J 


(A.  4) 


-  £  l£hi«i  -  ii’ 


plus  Eq.  (II. 10),  which  defines  d_^ ,  and  Eq.  (11.21), 


^  n  m  V  *  0* 
i-l  niH 


(11.21) 


We  discuss  k^,  the  thermal  diffusion  ratio  of  species  i,  in  Section  VI. B. 

Eq.  (A. 4)  makes  use  of  the  fact  that  each  multicomponent  diffusion  coefficient 
in  the  formulation  of  CC,  is  approximately  equal  to  the  first 

approximation  to  the  binary  diffusion  coefficient  .  We  have  discussed 


this  in  Section  III. 


The  equation  for  conservation  of  momentum  in  the  absence  of  external 


forces  is 


fr(pv  )  =  -  v  •  (pv  v  )  -  v  •  p. 
at  o  —  — o-o  — 

In  Eq.  (A. 5),  the  pressure  tensor  is 

l  •  PE  +  1, 


(A.  5) 


(A.  6) 


i 

where  IJ  is  the  unit  tensor  and  ! 


1  -  -  n  .  [Vv  +  (Vv  )T]  +  (2  .  —  k  ) (V  •  v  ) U . 

=*=  mix — o  ' — o'  1  v3  mix  mix  —  —o- 


(A.7) 


T 

In  Eq.  (A.  7),  nm^x  is  the  coefficient  of  viscosity  of  the  mixture,  (  )  is 
the  transpose  operation,  and  is  the  "bulk  viscosity"  coefficient.  The 

bulk  viscosity  is  closely  related  to  relaxation  effects  between  transla¬ 
tional  motion  and  the  internal  degrees  of  freedom  of  the  molecules  and  is 
equal  to  zero  for  mixtures  of  monatomic  gases.  Williams^  indicates  that 
Kmix  is  usually  negligible  for  combustion  processes. 

The  equation  for  conservation  of  energy  is 


M 

3t 


I  *  (Ev  )  -  V  •  (v  •  P)  -  7  •  a. 


(A.  8) 


where  E  is  the  energy  density  and  £  is  the  heat  flux.  The  equation  defining 
energy  density  is 
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in  which  e  is  the  internal  energy  per  unit  volume.  The  heat  flux,  £, 

^  ^  u  30,79 

is  given  by 


q  -  -  X°  .  VT  +  2  2  7Z.  7* 

mix—  .  ,  i  — 

i=l  ct 


h“  +  — 

i  n 


v 

Lt 


-vHCB 

nl^Ti  —  — 

frsr-,  (1,-vj. 


(A.  10) 


In  Eq.  (A.  10),  is  the  thermal  conductivity  of  the  gas  mixture  when  we 

ignore  the  transfer  of  energy  between  translational  and  internal  degrees 
of  freedom  of  the  molecules  and  the  effects  of  chemical  reactions.  The 
second  term  includes  the  transport  of  energy  by  excited  quantum  states 
of  gas  i  and  also  accounts  for  the  transport  of  heat  resulting  from  changes 
in  chemical  composition  with  temperature.  This  term  is  thus  directly 
related  to  the  molecular  internal  degrees  of  freedom  and  to  chemical  reactions. 
We  have  suppressed  the  index  i  of  for  brevity  and  have  used  °fK  to  denote 
total  number  of  moles  in  the  mixture  per  unit  volume.  We  also  point  out  that 
Eq.  (A. 10)  treats  each  quantum  state  as  a  separate  species;  thus  is  the 
diffusion  velocity  of  component  i  in  state  and  is  the  enthalpy  per  mole 
of  gas  i  in  the  quantum  state  a^.  The  last  term  in  Eq.  (A. 10)  represents  the 
"Defour  effect"  and  may  be  restated  in  terms  of  the  diffusion  ratios  {k.,^} 

(see  eq.  (A. 13)). 

Appendix  D  shows  how  we  may  separate  the  second  term  in  Eq.  (A. 10)  into 
two  parts — one  which  depends  on  the  quantum  states  {a^}  of  each  component  i 
and  the  other  which  depends  on  the  average  enthalpy  per  mole  of  each  component  i: 


%nd  Term 


int  ,  comp 

3  +  3 


v  v  _ 

=  S  E  £“  hJ  +  Z  <Hi>  . 

i*l  a  i*l 


(A. 11) 


The  quantity  is  the  molar  flux  of  component  i  in  quantum  state  ot^,  and 

is  the  enthalpy  per  mole  of  component  i,  averaged  over  the  quantum 

states  {a^}.  Thus  is  the  molar  enthalpy  which  we  would  expect  to 

obtain  experimentally  for  a  mixture  in  equilibrium.  Under  the  conditions 

l.nt 

listed  in  Section  V.A  of  this  paper,  of  E<j.  (A. 11)  becomes 

i 

ilnt  -2  SgJhJ.  -^mix  ^  (A*12) 

i«l  a 


where  then  represents  a  correction  to  the  thermal  conductivity  of  the 

gas  mixture.  We  derive  a  useful  equation  for  A*°^  in  Section  V.B. 

The  term  £Comp  in  Eq.  (A. 11)  represents  the  transport  of  heat  caused 
by  changes  in  chemical  composition  with  temperature,  and  under  some  condi¬ 
tions,  we  may  transform  £COmp  to  a  term  -A^°°P  VT.  Because  these  condi¬ 

tions  include  chemical  equilibrium  at  the  local  temperature  T,  which  does 
not  hold  in  general,  and  because  the  form  of  £COinp  in  gq.  (A. 11)  is  con¬ 
venient  for  calculations,  we  do  not  attempt  to  derive  A^°”p  in  this  paper. 
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By  comparing  terms  in  the  derivations  of  £  given  by  HCB  and  CC  , 
we  find  that  the  Defour  effect  may  be  represented  by 
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Appendix  B.  Sonine  Polynomials 

Here  ve  use  the  notation  of  HCB,  defining  the  Sonine 
polynomials  by  the  equation 


.(«>  , , .  £  ^ 

n  W  j.*o  (n+j  ) !  (m- j  ) !  j ! 


(B.l) 


with  the  orthogonality  condition: 


/*“  e‘*  S  <">(*)  s  d* 

*/0  n  n 


m:  mm 


(B.2) 


If  we  substitute  for  x,  and  use  n  »  3/2,  5/2,  eq.  (B.2)  gives  us 


f  f  (W  2)  V  2  d  V  *  — - —  5 

J  ri  &3/2  ^i  ;  i  a  H  m  mO 


(B.3) 


and 


I 


f.C0:i  si?2  (W,2)  V„4  d  Y±  «  15 n  i  (^)2  (B.4) 


5/2  v"i  '  vi 


(  J  5  (m) 

When  we  expand  A.  and  C.  in  terms  of  the  set  (S  }  and  B.  in  terms  of 
(ffl)  1  3/2  1 

an<i  then  substitute  these  expansions  into  the  integral  equations  for  the 

transport  coefficients,  we  obtain  expressions  of  the  form  of  eq.  (B.3)  and 

(B.4)  respectively.  For  this  reason,  the  Sonine  polynomials  constitute  a 

convenient  choice  for  calculating  approximate  values  for  the  transport 

coefficients. 
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Appendix  C.  First  Approximations  to  ^°mix 

3  7 

In  this  appendix  we  compare  the  values  of  derived  by  Mason, 

39  41 

Muckenfuss  and  Curtiss,  and  CC  ;  we  do  so  to  insure  that  we  can  derive 
eq.(V.12)  starting  with  any  of  these  "first  approximations"  and  to 
indicate  how  the  derivation  proceeds.  From  this  we  shall  see  that,  in 
eq. (V. 11) , 


CXmix(i)Jl  “ 


+  A 


v 

Z 

k-1 

Wi 


.  (C.l) 


*k  Gik 


The  importance  of  eq.  (C.l)  is  that  we  have  a  simple  method  of 
calculating  the  quantities  ,  which  we  need  in  order  to  compute 

approximate  values  of  the  thermal  diffusion  ratios  (Section  VI) . 


First  we  must  develop  equations  for  [ .  ],  from  the  formulation  of 

mix  i 

41  37 

CC  which  appear  in  the  same  form  as  those  of  Mason  and  Muckenfuss 

39 

and  Curtiss.  .  In  the  notation  of  CC,  we  have 


alS  tXraix(i)]l 


(C.2) 


which  are  solutions  to  the  set  of  coupled  linear  equations. 


Lii +  £  ht 


(1*1,2, «*«.,v). 


(C.3) 
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where 


L..  -  -± 


T*k 


*ii 


V 

+  2  5PfjD  ] 

lc-1  ^ikJl 

k*i 


6°1  +  s-4Bik>  mi  *  hiVn 


(mi +  \y 


L1J 


11-4B  .  -  8A. . 
_ ij _ il 

+  m^)2 


(C.4) 


And  A^j  and  are  ratios  of  collision  integrals  defined  by  CC.  Next 
we  solve  eq.(C.3)  by  applying  Cramer's  rule  to  obtain 


\l  Lv2 


column 

l 


'll 

L12  ‘  * 

•  •  •  •  • 

.  .  L. 

lv  | 

* 

* 

'21 

L22  *  * 

•  •  x^ •  •  • 

•  *  L2v 

vv 


(C. 


where 


[L  [  is  the  determinant  of  the  matrix  with  elements  {L^}.  We 


,th 


now  multiply  the  i  columns  of  the  numerator  and  the  denominator, 
respectively,  by  -4x^  for  each  i  *  l,2,...,v  and  factor  4  from  the  1 
column  of  the  numerator  to  obtain 


th 


column 

l 


a£  "  4 


L11 

^12  *  *  •  • ~XiX£ ’  • 

•  ,Liv 

L21 

L22  *  *  •  *~x2xr  * 

‘  *L2v 

• 

• 

• 

• 

i. 

• 

• 

• 

Lv2  *  *  *  *“W  * 

i  i*  1 

(C.6) 


where 


|l|  is  the  determinant  of  the  matrix  with  elements  (L^  given  by 


4xj  v  4TxiXk  (ll-4Bik)m2  +  gm^ 

i]l  ’  k>X  5P^ik]l  (miTV 

k*i 


Jii  “  [X 


ik 


(C.7) 


and 


4T  XiXi  mimi 

Lij  "  5P  TZrT'(n  )2  (11  -  4Bij  -  8Aij)- 
i*j  1  1  J 

By  eq.  (V.ll),  eq.  (C.2),  eq.  (C.6)  and  eq.  (C.7),  we  then  have 


ix\  ]CC 

mix  1  ■  4 


L11  L12 
L21  L22 


Lvl  Lv2 


X1  X2~ 


•Llv  X1 
*L2v  X2 


.  L  x 

VV  V 

.x  0 

V 


(C.8) 


where  the  notation  [  ]  cc  identifies  the  formulation  as  that  of  CC. 
Eq.  (C.8)  is  in  the  same  form  as  the  corresponding  equations  of 

’  4 

Mason  and  of  Muckenfuss  and  Curtiss  the  elements  {L^}  differ, 
however.  To  compare  the  respective  sets  {L^},  we  use  the  following 
relations  in  eq.  (C.7): 
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a. 


M, 


m.+m. 
i  J 


M.+M. 
i  3 


(i,j  -  1,  2 


v  ) 


CC.9) 


and 


84,  85 


I  a*HCB 
5  ij 


3ec 

ij 


■  —  B  *  HCS 
5  ij 


CCC-1953  m  i  +  CCC-1970  ,  6  C*HCB 


(C.10) 


The  quantity  C„  is  another  ratio  of  collision  integrals  which  we 
encounter  in  Section  VI,  and  the  notations  CC-1953  and  CC-197Q 
differentiate  between  the  second  and  third  editions  of  CC.  We  then 
have,  from  the  formulation  of  CC, 


‘ii 


4*i  _  16T  r  *i*k  {(  4  ~3Bik)Mi  *  ^i^ik1 

[X“]L  25 P  k=1  t^kl1 

k*i 


16T  xx,  MM,  ==  * 

t  m  _  1  J  ,  Ji  j  r jj.  -  3B  -  4A  \ 

Lij  25P y). ,  J .  (M  +M  )2l4  Jfiij  Aij^* 

i#j  1  1  J 


(C.ll) 


The  corresponding  equations  derived  by  Mason 
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are 


4*l  161  }  Vk  (TM1 -1-? 

11  -  -  1X-],  -  25P  (M±  *  Spr- 

k*i 


and  (C . 12 ) 

,  16T  _Vl_  JS _  {10  .  4a*  } 

Lij  25P  t0.,L  (M,+M.)2  uu  Aij;’ 
i+j  13  1  1  J 


and  the  matrix  elements  of  Muckenfuss  and  Curtiss' 


.14,  39 


are 
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ii_  is  ?  w 

•  '  25P  t£1  ^17  (Mi  +  V2 

k*i 


and 


,  16T  _Vj_  JVJi _  {55  _  * 

ij  25P  0  }  (M  +M  )2l4 

i/j  13  1  13 


4A*tj} 


(C.13) 


Notice  that  the  only  differences  among  the  three  sets  of  results  occur 

*  it 

in  the  bracketed  portions  concerning  the  integrals  A  and  B  ;  in  fact, 

L^Ci^j)  is  the  same  in  eq.  (C.ll)  and  eq.  (C.13).  To  derive  the 

approximate  expression  in  eq.  (V.12)  according  to  the  method  of  Mason 
14 

and  Saxena,  we  would  set  L^(i#j)  *  0  in  eq.  (C.  11-13).  We  would  then 
represent  the  portions  of  L,^  containing 

empirical  constant.  Since  the  remaining  portions  of  eq.  (C. 11-13) 
are  identical,  all  three  formulations  will  lead  to  the  same  semi- 
empirical  equation. 

Using  the  above  procedure  on  eq.  (C.5),  we  can  see  that  eq.  (C.2) 

gives 


V 


Aik’  and  Btk  by  “ 


[Xmix(i)3l  "  *  (C. 14) 

14 

Continuing  with  the  treatment  of  Mason  and  Saxena,  the  reader  can 
confirm  that  eq.  (C.l)  is  correct.  As  we  indicated  previously,  this  is 
important  in  calculating  the  coefficient  of  thermal  diffusion,  since  we 
have  a  simple  method  of  evaluating  which  appears  in  the 

approximate  expression,  eq.  (VI.14). 
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Appendix  D.  Decoupling  the  Effects  of  Molecular  Internal  Degrees  of 
Freedom  and  Chemical  Reactions  in  the  Heat  Flux  Equation 

In  Hirschf elder's  treatment  of  the  thermal  conductivity  of  a  multi- 
30 

component  gas  mixture  (see  Section  V),  he  has  discussed  the  heat  fluxes 
<^nt  and  the  equations  for  which  depend  on  the  molecular  internal 

degrees  of  freedom  and  the  chemical  composition  of  the  mixture,  respectively. 
In  doing  so,  however,  he  did  not  show  how  to  decouple  the  terms  representing 
these  effects  in  eq.  (A. 10)  for  the  total  heat  flux  £.  Instead  he  assumed 
that  eq.  (A. 11)  had  already  been  derived  and,  therefore,  that  he  could 
consider  each  factor  independently.  Below  we  show  how  to  derive  eq.  (A. 11) 
by  defining  the  heat  fluxes  £  Qt  and  ^comp  in  a  manner  consistent  with 
Hirschf elder 's  discussion.  In  Appendix  E,  we  demonstrate  that  our  expres¬ 
sions  for  £int  and  ^comP  in  fact  decouple  heat  transport  caused  by  changes 
in  the  population  of  quantum  states  and  the  chemical  composition  with  tem¬ 
perature  ,  and  we  show  that  Hirschf elder ' s  discussions  are  rigorously  valid 
for  chemically  reactive  mixtures  of  polyatomic  moldcules. 

Our  derivation  begins  with  Eq.  (A-10)  for  the  total  heat  flux  g^,  which 
gives  us 


a 


-  X°.  VT 
mix  — 


v 

+  2 

i«l 


(D.l) 


in  the  absence  of  thermal  diffusion.  In  Eq.  (D.l), 


a  yO 
i  ~i 


(D.2) 
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is  the  number  of  moles  of  component  i  in  quantum  state  a^  which  are  moving 
in  1  sec  through  a  1  cm2  surface  which  is  perpendicular  to  the  velocity 
V?»  and  is  the  enthalpy  per  mole  of  component  i  in  quantum  state 
The  quantity  is  the  diffusion  velocity  of  component  i  in  quantum  state 
given  by 


V.“*  v,a-  v  ,  (D. 3) 

— i  —  i  — o 

and  x^1  is  the  mole  fraction  of  component  i  in  quantum  state  a^, 

! 

a 
n. 

a  i 
x.  *  — 
i  n 

with 

a  3  S  S  n“ 
i-1  a 

We  have  left  the  index  i  off  a^  in  all  equations  for  brevity.  Notice  also 
that  these  and  many  subsequent  equations  are  restatements  of  those  in  Section 
II  assuming  that  each  combination  (i,  a^)  is  a  separate  component  of  the 
mixture.  Thus  we  have  for  the  fluid  velocity  v^  the  relation 

S  S  (D.6) 

i-1  a 

which  gives  us  the  equation 


v 

—o 


s  ; 

i-l 


a 

n  m 
li 


given  by 

(D.4) 

(D.5) 
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(D.  12) 


and 


comp  _  V-. 

£  -  L 

i-1 


V 

i-J 

a 


(D.  13) 


Next  we  define  the  molar  flux 


(D.14) 


where  is  given  by 

I 

-  V*  -  ±  Snf  m  ,  (D.15) 

^  *  P±  g  1  1 

the  average  velocity  of  component  i  in  quantum  state  ot  relative  to  the  mass 
average  velocity  of  component  i.  From  Eq.  (II. 6)  and  Eq.  (D.15),  we  obtain 


where 


(D.16) 


<?i>=^  ~  “i* 


$  -  6 


(D.17) 


with  the  notation  <  >  signifying  an  average  over  the  quantum  states  of 
species  i.  We  note  that,  in  equations  (D.15)  and  (D.17),  we  have  used 


v 


(D.18) 


and 
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(D.19) 


pi  "  ni  °i  "  ?  ni  V 

P 


From  Eq.  (D.14)-(D.18) ,  we  obtain 


S  G“  =  0,  (D.20) 

— i 
a 

which  is  the  condition  used  by  Hirschfelder^  to  derive  for  a 

mix 

nonreacting  gas  mixture  in  a  steady  state.  Thus  Eq.  (D.ll)  will  give  the 
same  formula  for  in  the  more  general  case  of  a  reactive  gas  mixture. 

Now  we  must  show  that  our  definition  of  £COmp  in  (D.12)  is  compatible 
with  Hirschf elder's  derivation  of  the  heat  flux  resulting  from  changes  in 
chemical  composition  with  temperature.  From  Eq.  (D.2),  (D.3),  (D.13),  and 
(D.16),  we  have 


ii*  m^xi[3  -  -  <5  -  <v  *  v,) 

From  Eq.  (D.6  )  and  (D.17),  we  see  that 

D  n  ®,(<v.>-  v )  -  0.  (D.22) 

i-1  1  1  ^  ^ 

Further  the  total  molar  flux  of  species  i  in  a  reference  frame  moving  at 
the  velocity  is  given  by 

-S  g“  "  2  (fij  +  "  S  ^°-Hx  (<  vi>  -  yo)  .  (D.23) 

a  a  a 
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From  Eq.  (11.20)  and  the  accompanying  discussion  of  diffusion  velocities 


,  we  also  have 


*i  m7l*  i  ^i 

(D.24) 

Thus  from  Eq.  (D.22)-(D. 24) ,  we  obtain 

(D.25) 

and 

,<x  an  a  ~  ! 

%  ’7*1  4' 

(D.26) 

Using  Eq.  (D.26)  in  Eq.  (D.12),  we  find  that 


lco"p  -  S  V,  n“ 

i-1  a 


(D.27) 


*  <“i>  • 

where  as  in  Eq.  (D.17),  we  have  defined  (  H^)  as  the  enthalpy  per  mole 
averaged  over  the  quantum  states  a^. 


<V  "r*  Sxi  K- 

3  *i  a  1  1 


(D.28) 


We  interpret  <U^)  as  the  enthalpy  which  is  measured  experimentally;  thus 
Hirschf elder *s  assumption  that  the  gas  is  in  equilibrium  at  the  local 


temperature  T  is  necessary  if  we  use  tabulated  values  of  the  enthalpy  in 
our  models. 

30 

Hirschfelder  used  an  equation  identical  in  form  to  Eq.  (D.27)  to 
discuss  the  effect  of  changes  in  chemical  composition  with  temperature  upon  the 
thermal  conductivity  of  a  gas  mixture;  again,  however,  his  treatment  contains  the 
implicit  assumption  that  the  terms  associated  with  molecular  internal 
degrees  of  freedom  will  not  affect  the  results.  In  Appendix  E,  we  outline 
Hirschfelder 's  discussion  and  using  the  equations  which  we  have  developed 
in  this  appendix,  we  show  that  the  above  assumption  is  correct. 
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Appendix  E.  Derivation  of  A and  A0*?®** 

mix  mix 

Below  we  outline  Hirschf elder *s  derivation  of  A*®fc  and  AC°m^  from 

mix  mix 

<l^nt  and  ^comP >  respectively.  In  our  discussion,  we  will  use  the  more 
general  framework  of  Appendix  D  to  demonstrate  the  validity  of  his  treat¬ 
ment.  As  an  example,  we  will  then  investigate  in  detail  the  case  of  a  gas 
with  only  one  chemical  species. 

30 

To  simplify  our  discussion  and  to  parallel  Hirschf elder's  derivation, 
we  will  assume  that  a  temperature  gradient  exists  only  along  the  z-axis,  that 
no  heterogeneous  reactions  take  place  (i.e.,  that  the  bounding  surfaces  are 
chemically  inert  and  have  no  catalytic  activity) ,  and  that  the  system  has 
reached  a  steady  state.  The  last  assumption  means  that 


(E.l) 


for  each  species  i,  and  eq.  (A. 2)  becomes 


fe  <wv> 


«i- 


W 


(E.2) 


Following  Hirschf elder,  we  transform  to  a  local  reference  frame  in  which 
the  fluid  velocity  vq  is  zero.  We  may  then  express  eq •  (E.2)  in  the  form 


3z 


V 


(E.3) 
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where  is  the  net  number  of  moles  of  component  i  formed  per  cm  per 
second.  Substituting  eq.  (D.23)  into  eq.  (E.3),  we  have 


^®i  3  r-v  a  3  v'  -0* 

3z~  =  Iz  ^  gi  *  3z  ^  gi 
a  ai 


V 


(E.4) 


and  we  see  that  only  the  part  of  the  molar  flux  which  is  associated  with 
^comp  ^  eq.  (Q^2)  is  related  to  chemical  reactions,  as  we  would  expect. 
To  transform  eq.  (D.ll)  and  (D.12)  to  expressions  of  the  form 


a  -  -  *£t. 


(E.5) 


we  assume  that  thermal  and  pressure  diffusion,  radiation  processes,  and 
external  forces  are  negligible.  Under  these  conditions,  eq.  (A. 4)  becomes 


2 

dz 


S  2 

j«l  8 


a  8  6  a 

Xi  g1  -  X1  8i 


«n«oa0 

ij 


(E.  6) 


in  which  we  distinguish  between  molecules  in  different  quantum  states. 

We  have  suppressed  the  notation  [  ] ^  in  representing  the  first  approxi- 

mation  toJO^J,  the  binary  diffusion  coefficient  for  component  i  in  state 
and  component  j  in  state  8^.  We  may  now  separate  the  effects  of  internal 
degrees  of  freedom  from  those  of  chemical  reactions  by  substituting  Eq.  (D.13) 
into  Eq.  (E.6)  to  obtain 
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where 


3 

dz 


dx" 


'<] 

int 

.  a* 
dxi 

dz 

+ 

dz 

t- 

•»  m 

-»  comp 


dz 


V 

=  s 

J-l  8 


xi  G?  -  G? 
V  i  1  1  i 


*>\t& 


and 


dx. 


comp 


v 

=  s  s 

J-l  8 


liL 


ctB 

ij 

-  z'1 

X1  8i 


(E.7) 


(E.8) 


o8 

ij 


(E.9) 


In  Eq.  (E.8)  and  (E.9),  we  have  again  suppressed  the  dependence  of  and  8. 
on  i  and  j,  respectively.  Sunming  Eq.  (E.8)  over  cx^,  and  referring  to 
Eq.  (D.20),  we  see  that 


s 

a 


0 


(E.10) 


so  that  the  gradient  in  chemical  composition  is  given  by 


dx^ 

dz 


r  0n  comP 

dx. 


2 

a 


dz 


(E. 11) 


Eq.  (E.ll)  demonstrates  the  natural  decoupling  of  the  effects  of  internal 
degrees  of  freedom  and  chemical  reactions  which  the  framework  of  Appendix  D 
permits.  Performing  the  sum  over  a  in  Eq.  (E.ll)  and  noting  Eq.  (D.23),  we 
find  that 


Vi’ 


x1*i 


^j 


(i  -  1,  2,  •••,  v). 


(E.12) 
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where  we  have  dropped  the  dependence  of  the  binary  diffusion  coefficients 

30 

on  quantum  states  a^and  Sj following  Hirschf elder 1 s  assumption.  (We  have 
discussed  the  assumption  that  is  independent  of  and  8^.  in  Section  V.A.) 

If  we  further  assume  chemical  equilibrium  at  the  local  temperature  T, 

Gq.  (G.3)  becomes 


Ri^Xl’  X2’  V  T* 


(E.13) 


from  which  we  may  obtain  a  set  (x^e(T,P)}  describing  the  equilibrium  compo¬ 
sition  of  the  system.  Hirschfelder"^  indicates  that  these  will  provide  a 

first  approximation  to  the  actual  composition.  We  may  then  solve  Eq.  (E.12) 

30 

for  the  set  {g^  to  obtain 


_  [«i  ~  „HCB  field! 

L » >i 1J  J  di  J di 


(E.14) 


where  is  a  multicomponent  diffusion  coefficient  (formulation  of  HCB) 

and  the  subscript  "e"  indicates  that  we  have  substituted  the  set  (xie(T,P)} 
into  our  solution.  From  Eq.  (D.24)  and  (D.27),  we  see  that 


^comp  dT 
'mix  dz 


where 


5  v  v  ax. 

comp  Sv  dHCB  <H  , 

mix  0  i-lj-1  iJ  J  1  dT 


(E.15) 


(E. 16) 
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For  further  discussion  of  Eq.  (E.16),  we  refer  the  reader  to  Hirschfelder 

80 

and  Butler  and  Brokaw  . 

We  now  return  to  Eq.  (E.8)  to  derive  an  equation  analogous  to  Eq.  (E.5) 
and  (E.15)  for  the  thermal  conductivity  related  to  internal  degrees  of  freedom. 
Performing  the  sum  over  B  and  using  Eq.  (D.4),  (D.18),  and  (D.20),  we  obtain 


p  -i  int  .  . 

dx“  x.  v  x.  \ 


(E.16) 


in  which  we  have  again  assumed  that  the  binary  diffusion  coefficients  are 
approximately  the  same  for  all  and  8^ .  If  each  component  i  is  in 
equilibrium  at  the  local  temperature  T,  the  mole  fractions  x”  are  functions 
of  local  temperature  and  pressure,  and  we  have 


(E.17) 


We  may  now  use  Eq.  (E.16)  and  (E.17)  in  Eq.  (D.ll)  to  obtain 


in  which 


int  .int  dl 
q  *  ’  \nix  dl 


(E.18) 


(E. 19) 
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To  complete  the  discussion,  we  will  now  compute  A^nt  for  the  case  of 
a  gas  consisting  only  of  component  i.  From  Eq.  (E.19),  and  the  fact  that 


P  l‘J 

2  (  0  i  t  j 


(E. 20) 


we  have 


,  int 


dx 


VMi-  dT“i 

a 


(E. 21) 


By  Eq.  (D.28),  we  know  that 


<H.  ) 


a 

i-j  xi  ’ 


(E. 22) 


and  the  molar  specific  heat  at  constant  pressure  for  the  gas  i  is 


d  <H±  ) 


.a 


a 


pi  dT 


d  dx, 

id!  +  L  Hi  F 
a  a 


=  ^  xa 


(E.23) 


The  first  term  is  the  "translational"  contribution  to  the  heat  capacity 
30 

and  equals  5p/2  .  Thus  we  find  that 


-  !  *> 


(E. 24) 
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j 


and  the  total  thermal  conductivity  for  species  i  becomes 


Ai  = 


a?  +  . 

i  wii  pi 


(E.25) 


Defining  6  by  the  equation 


5 


fi 


5R  . 
11 


(E.26) 


we  obtain 


\°  *  1  "  6fi  +  5  6fi  “R~  =  Ei^fi^  * 


(E.27) 


Using  detailed  calculations  with  the  Buckingham  and  Lennard-Jones  potentials, 
29 

Hirschfelder  has  calculated 


Sfi  *  0.885. 


(E.28) 


This  gives  us  the  value  of  Ei  shown  in  Eq.  (V.8). 
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HCB  CC 

Appendix  F.  Demonstration  that  kTi  ■  kTi  for  Binary  Mixtures 

In  this  appendix,  we  use  the  formulation  of  kinetic  theory  given 
CC 

by  CC  to  calculate  for  a  binary  gas  mixture,  and  we  compare  the 

HCB 

result  to  the  expression  for  k^  .  We  do  this  to  demonstrate  the 

following:  that  (a)  different  formulations  will  give  values  of  k, 

CC 

which  are  equal,  (b)  A^.  and  D^.  are  not  equal,  and  (c)  eq.  (VI. 6) 

CC  HCB  W 

relating  D^,  DTi  ,  and  is  valid.  We  calculate  A  in  terms  of 

JS  «•  which  is  equal  to  for  binary  gas  mixtures,  from  the  equation 


Ti 
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where 


V 

Y 

'l2 

V2 
JO  12 

*1*2 

V 

$  21 

JO  22 

P1 

p2 

(F.l) 


(F.2) 
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and  Y  is  the  cofactor  of  ~  ^  in  Q. 
13  X>ij 

of and^^,  we  °*5tain 


From  eq.  (F.l-2)  and  the  equality 


(F.3) 


»  -ex  +  »j. 


and  the  above  equations,  we  find  that 


_CC  °2  +  °lp2  .CC 

Dn  •  rz -  Si 


,CC  „ 
12  kTl  P2 

X1X2P 


We  note  that  eq.  (II. 6)  and  eq.  (VI. 6)  give  us 


Px  •  n1m] ,  P2  -  n2m2’ 


(F.8) 


(F.9) 


(II. 6) 
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and 


„  nCC  HCB 
P1  °T1  "  DT1  * 


Substituting  these  equations  into  eq.  (F.9),  we  obtain 


.CC  , 

.HCB  w  12  kTl  °  ala2 
T1  "  o 


CC 


Solving  for  k^  and  comparing  this  with  the  corresponding 
for  k!^®,  we  find  that 


„HCB 

CC  _  °T1  p 


HCB 


^Ti  n2  ®L®2  kn  • 


(VI. 6) 


(F. 10) 
equation86 


(F.ll) 
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